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Many physical systems can be modeled as large sets of domains “glued” together along bound¬ 
aries - biological cells meet along cell membranes, soap bubbles meet along thin films, countries 
meet along geopolitical boundaries, and metallic crystals meet along grain interfaces. Each class of 
microstructures results from a complex interplay of initial conditions and particular evolutionary 
dynamics. The statistical steady-state microstructure resulting from isotropic grain growth of a 
polycrystalline material is canonical in that it is the simplest example of a cellular microstructure 
resulting from a gradient flow of a simple energy, directly proportional to the total length or area 
of all cell boundaries. As many properties of polycrystalline materials depend on their underlying 
microstructure, a more complete understanding of the grain growth steady-state can provide insight 
into the physics of a broad range of everyday materials. In this paper we report geometric and topo¬ 
logical features of these canonical two- and three-dimensional steady-state microstructures obtained 
through large, accurate simulations of isotropic grain growth. 

PACS numbers: 61.72.-y, 61.43.Bn 


I. INTRODUCTION 

A polycrystalline material is an ensemble of individ¬ 
ual crystallites or grains, joined along grain boundaries. 
Polycrystalline microstructures are often formed during 
material synthesis and continue to evolve through pro¬ 
cesses such as grain growth and recrystallization. Most 
metals and ceramics, and many semiconductors, ubiq¬ 
uitous in every day use and technology, are polycrys¬ 
talline. One of the most common class of polycrys¬ 
talline microstructures is that resulting from normal 
grain growth. Normal grain growth is the process in 
which grain boundaries migrate with a velocity propor¬ 
tional to their mean curvature, leading to a monotonic 
decrease of grain boundary area and excess energy, and a 
monotonic increase in mean grain size. In its ideal form, 
the grain boundary energy and mobility are isotropic 
and constant. It has been observed that normal grain 
growth leads to the evolution of microstructures that ex¬ 
hibit a statistical steady-state, during which geometrical 
and topological properties of the microstructure remain 
statistically constant as the mean grain size increased. In 
many real experimental cases this picture is complicated 
by anisotropy and the presence of solutes or additional 
phases. Nonetheless, isotropic normal grain growth mi¬ 
crostructures are the most generic of all polycrystalline 
microstructures. Because of the simplicity of the grain 
growth process in the isotropic limit, and the ubiquity of 
grain growth microstructures it produces, normal grain 
growth microstructures can be thought of as the “har¬ 
monic oscillators” of polycrystalline microstructures or 
“the canonical polycrystalline microstructure.” 

We note that other classes of generic polycrystalline 
microstructures exist. Eor example, random grain nucle- 
ation and constant grain boundary velocities, appropri¬ 


ate for modeling certain solidification and phase trans¬ 
formation processes, leads to microstructures which are 
Voronoi tessellations of Poisson point processed. These 
Voronoi microstructures have fiat grain boundaries, and 
distributions of their geometric and topological features 
are vastly different from those commonly observed in real 
polycrystalline materials; we return to this comparison at 
the end of this report. 

Characterizing steady-state microstructures experi¬ 
mentally in three dimensions is challenging due to the 
inherent difficulty in exploring the geometry and topol¬ 
ogy of grains buried deep within a three-dimensional mi¬ 
crostructure. Several such studies have been performed, 
including by embrittling the polycrystal grain bound¬ 
aries and subsecmently separating individual grain^, by 
serial sectioning^^^, and more recently by synchrotron- 
based diffraction technique^^^^^. Because of the difficul¬ 
ties in experimentally characterizing three-dimensional 
microstructures consisting a large numbers of grains, and 
the fact that isotropic grain boundary energy and mobil¬ 
ity are rare, simulations provide the preferred method 
to investigate isotropic normal grain growth microstruc¬ 
tures. 

The primary purpose of this paper is to provide 
the most complete characterization of the geometric 
and topological features of two- and three-dimensional 
isotropic grain growth microstructures to date, as well 
as the correlations between such features. To th is end, 
we employ a computer simulation techniqu d^^ l ^^ l to gen¬ 
erate one of the largest and most accurate steady-state 
isotropic grain growth microstructure databases. Our ge¬ 
ometric and topological characterization of these ideal, 
normal grain growth microstructures provides a rigorous, 
fundamental set of data against which all polycrystalline 
microstructures and grain growth theories can be com- 
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FIG. 1. (a) A two-dimensional grain-growth microstructure, (b) a cross-section from a three-dimensional grain-growth mi¬ 
crostructure, and (c) a three-dimensional grain-growth microstructure. 


pared. 

A secondary focus of this paper is to compare two- 
and three-dimensional steady-state microstructures in or¬ 
der to explore the manner in which the introduction of 
a third dimension essentially changes the steady-state 
structure. This is partly motivated by the substantial in¬ 
vestment of the research community in two-dimensional 
simulation^^^^^, often with the implicit hope that re¬ 
sults gleaned from these will shed light on the experimen¬ 
tally more relevant three-dimensional case. Finally, since 
three-dimensional experimental samples are often viewed 
and characterized in cross-section, we also examine the 
relationship between cross-sections of three-dimensional 
structures with both the two- and full three-dimensional 
microstructures. 


II. SIMULATION 
A. Prior approaches 

Prior simulations of grain growth fall rou ghly into 
three categories. Monte Carlo Potts models^^EHIlIini 
were among the earliest and are still among the most fre¬ 
quently implemented. Here, space is subdivided into cu¬ 
bical voxels which are assigned a label indicating to which 
grain they belong. Voxels are relabeled in a stochastic 
manner to minimize interfacial energy between neighbor¬ 
ing grains. While the implementation of this method is 
relatively straightforward, some of the compli cation s that 
arise include anisotropic boundary condition^^^^, weak¬ 
ening of vertex angle bounda ry condi tion^^^^, boundary 
pinning, and altered kinetics^^I^^^. Closely related to 
this model is the cellular automaton modeP^^^, which 
allows additional flexibility in defining energetics. 

Phase-field mod els co nstitute a second category of 
simulation methodJ^^^^. In this approach, grains are 
represented as regions in which the value of a nom¬ 
inally continuous order parameter goes to one. The 
evolution of the microstructure is described using time- 
dependent Ginzburg-Landau equationJ^, where differ¬ 


ent order parameters are associated with each grain ori¬ 
entation. Modifications of the energy functional allows 
for the incorporation of additional physical effectJ^. 
Limitations of this approach include numerical instabil¬ 
ity and the high cost o f of representing systems with 
large numbers of grain^^^^. Variat ions on this ap¬ 
proach include continuum-fielcP^^^^^^ and multi-phase- 
fielcP^ models. The simulation of three-dimens ional 
grain growth by diffusion-controlled interface motiorP^^ 
shares some features with phase-field models, and ap¬ 
pears to be unconditionally stable and conceptually and 
computationally straightforward, allowing for the simu¬ 
lation of a large number of grains. 

A final category of grain growth simula tions are front¬ 
tracking methods, including vertex model4^^^, finite el¬ 
ement models^^M^ and d iscret ized boundary models as 
developed by the author^^SHIl employed here. In 
such front tracking approaches, the grain boundary net¬ 
work is explicitly discretized and evolved without devot¬ 
ing memory or processing power to the grain interiors. 
While computationally efficient and convenient for the 
measurement of geometric and topological features, this 
approach requires that all topological events which occur 
during normal grain growth be anticipated and handled 
explicitly. The complexity associated with mesh man¬ 
agement in three dimensions has been a limiting factor 
in the widespread implementation of these models. 


B. Current approach 

Most front-tracking algorithms evolve a system in the 
direction of steepest descent of an energy proportional to 
the total grain boundary aresP^l^, which is equivalent to 
mean curvature flow in a conti nuous setting. One notable 
feature of our simulation j^^ l ^^ Hs that the equations o f mo- 
tion derive directly from the von Neumann-MullinP^^^ 
relation in two dimensions and the MacPherson-Srolovitz 
relatioiPIl in three dimensions. This means that our sim¬ 
ulations satisfy the constraints imposed by these exact 
relations to high accuracy. The low curvature of most 
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grain boundaries justifies our use of a restricted dis¬ 
cretization to reduce the computational requirements. 

In two dimensions, boundaries that separate two neigh¬ 
boring grains are referred to as edges^ and three edges 
meet at a vertex. Each edge is discretized into a piece- 
wise linear curve; the discretization adapts during the 
microstructure evolution to ensure stability and numeri¬ 
cal accuracy. Equations of mo tion a re obtained from the 
von Neumann-Mullins relatioiSESl, which describes the 
evolution of individual grains: 

(>-?)■ « 

where A is the area of a grain, n is its number of edges, 
and M and 7 are scalar constants describing the grain 
boundary mobility and energy, respectively. This result 
is exact for isotropic grain growth in two dimensions, 
where the total energy of the system is proportional to 
the sum of all edge lengths and the system evolves in the 
direction of steepest descent in energy. This equation 
may be adapted to give equations of motion for a dis¬ 
cretized representation of the systenP^, which is evolved 
using an explicit time integration scheme. Since this al¬ 
gorithm derives from the exact von Neumann-Mullins 
relation, the behavior of every grain in the system sat¬ 
isfies this condition with minimal error. Specifically, the 
error resulting from the use of a discrete time step At is 
of order O(At^), and the error from the discretization of 
grain boundaries is small as welW Occasionally, when 
an element of the system becomes sufficiently small, we 
adjust the network topology by removing an edge or a 
grain, or by flipping an edge. Interested readers are re¬ 
ferred tcP^ for a more detailed discussion of the method 
and error analysis. 

The situation is similar in three dimensions, where a 
grain boundary separating two neighboring grains is re¬ 
ferred to as a faee. Three adjacent faces meet along an 
edge and four edges meet at a vertex. In simulations, the 
grain boundary network is discretized by triangulating 
grain faces, and edges are discretized into piecewise lin¬ 
ear curves. The extension of the von Neumann-Mullins 
relation to three dimensions is the MacPherson-Srolovitz 
relatiorl^. The relation describes the volume evolution 
of individual grains as 

where V is the volume of grain, £ is a one-dimensional 
measure of the grain called the mean widths and M. is 
the sum of lengths of all grain edges. This result is an 
exact description of isotropic grain growth in three di¬ 
mensions. As before, this equation may be adapted to 
give equations of motion for a discretized representation 
of the systenP^. An explicit time integration scheme is 
used to evolve the system. Since the equations of mo¬ 
tion are derived from the exact MacPherson-Srolovitz 
relation, the evolution of each grain volume satisfies this 


condition with minimal numerical erroi^m. Occasionally 
the size of a grain, face, or edge will become sufficiently 
small, requiring that the topology of the system be ad¬ 
justed appropriately. Interested readers are referred tcl^ 
ancP^ for a more detailed analysis of the method. 

One difficulty in using simulations to characterize the 
steady-state is that the system must be evolved from the 
initial system to the point that steady-state is achieved. 
This requires that the initial model contain several times 
more grains than eventually contribute to the measured 
steady-state microstructure statistics. Moreover, since 
the steady-state condition is not precisely defined, iden¬ 
tifying the point in time when it is reached is still a mat¬ 
ter of some contention. Data for our 2 D system came 
from a simulation which began from a Voronoi tessel¬ 
lation of 10 , 000,000 points uniformly distributed in the 
unit square with periodic boundaries. Data for our 3D 
system came from 25 simulations, each of which be¬ 
gan from a Voronoi tessellation of 100,000 points uni¬ 
formly distributed in the unit cube with periodic bound¬ 
aries. This number of grains is well beyond the range 
of contemporary front-tracking models^, and rivals the 
largest grai n growt h simulations performed to date by 
any methocP^^^^El. 


C. The Canonical Steady-State 

One reason that steady-state grain-growth microstruc¬ 
tures are so important is that despite the fact that the 
evolution equation (mean curvature flow) is very sim¬ 
ple, the resultant microstructures are particularly robust; 
i.e., a wide range of initial microstructures all eventually 
evolve to a statistically identical microstructure. To il¬ 
lustrate this point, we constructed three very different 
initial microstructures and allowed them to evolve for a 
time sufficient for their statistical properties to stabilize 
(i.e., reach the steady-state behavior). 

Each of the three initial conditions considered here 
were constructed as Voronoi tessellations of a set of points 
in space. The first system of points we considered re¬ 
sulted from a random Poisson process with 100,000 points 
in the unit cube; i.e., the x, and 2 : coordinates of each 
point were chosen with uniform probability in the domain 
[0,1]. To construct the second system, we generated a set 
of points on a 30 X 30 X 30 simple cubic (SC) lattice in the 
unit cube, for a total of 27,000 points. Then, each point 
was given an independent, random displacement from its 
lattice position according to a three-dimensional Gaus¬ 
sian distribution with a standard deviation equal to one 
tenth the SC lattice parameter. The third system was 
constructed in a similar manner as for in the perturbed 
SC case, except that the original coordinates were chosen 
to he on a 20 X 20 X 20 face-centered cubic (ECC) lattice 
with 4 points per unit cell, for a total of 32,000 points; 
in this case the points were perturbed away from their 
ECC positions according to a three-dimensional Gaus¬ 
sian distribution with a standard deviation equal to one 
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FIG. 2. (a) The evolution of the average faces per grain (F) in the three systems (circles indicate (F) at t = 0), converging to 
the steady-state value 13.766. (b) The evolution of the fraction of faces with 4, 5 and 6 edges in the three systems, converging 
to steady-state values p(4) = 26.2%, p(5) = 36.4%, and p{6) = 22.0%. (c) Initial (squares) and final (circles) distributions of 
grain volumes (normalized by the mean volume) in the three systems. 



FIG. 3. Gross-sections of three distinct initial microstructures 
(top) and the steady-state grain growth microstructure result¬ 
ing from it (bottom). Each system began as a Voronoi tes¬ 
sellation from a set of points: (a) random distributed points, 
(b) a perturbed simple cubic lattice, and (c) a perturbed face 
centered cubic lattice. 


fifth of the FCC lattice parameter. In all three cases, 
we evolved the systems according to the above algorithm 
for a sufficiently long time (measured in dimensionless 
units of I/M 7 ) that the statistical properties all achieved 
steady-state to within statistical error. 

Figure [^a) shows how the average number of faces 
per grain (F) changes as the three systems coarsen un¬ 
der curvature flow. The initial value of (F) is differ¬ 
ent in the three systems; in particular {F)py = 15.535, 
{F)sc = 15.693, and {F)fcc = 14.654 (as indicated by 
the symbols at t = 0). Over time, these values gradually 
decrease and all converge to 13.766 faces per grain. In a 
similar manner, Fig.j^b) shows the evolution of the frac¬ 
tion of faces with n = 4, 5 and 6 edges p{n) in the three 
systems. Here again we see that despite significant dif¬ 
ferences in their initial values, p{n) for all three systems 
converge to the same values, p(4) = 26.2%, p(5) = 36.4%, 


and p{6) = 22.0%. Finally, Fig. |^c) shows the initial 
and steady-state grain volume distributions, normalized 
by the mean grain volume, for all three systems. De¬ 
spite having started with drastically different grain size 
distributions, all three systems evolve to the same steady- 
state. The statistical properties of the three systems are 
effectively identical after t = 0 . 002 , supporting the no¬ 
tion that initial conditions are eventually “forgotten” af¬ 
ter a system has coarsened by grain growth for a sufficient 
time. 

The convergence of different random microstructures 
to the same statistical steady state can also be observed 
directly; typical initial and final cross sections from the 
three types of initial conditions are shown in Figure. 
Notice that despite the different initial appearances and 
degrees of order, the final three microstructures could 
easily be different areas of a single structure. 


III. CHARACTERIZATION 


This section reports on data from microstructures ini¬ 
tialized as Poisson-Voronoi (PV) cells (i.e., Voronoi tes¬ 
sellations of randomly distributed points) and evolved by 
the grain growth (curvature flow) algorithm until a statis¬ 
tical steady-state was achieved. A single two-dimensional 
simulation began with 10 , 000,000 grains and achieved 
steady-state when approximately 4% of the grains re¬ 
mained. Samples taken from several points in the subse¬ 
quent evolution gave two-dimensional statistics from a to¬ 
tal of 1,000,000 grains. In three dimensions a series of 25 
simulations was performed, each beginning with 100,000 
grains. Statistical steady-state was achieved when just 
over 10 % of the grains remained, and samples taken 
from the resulting configurations gave three-dimensional 
statistics from a total of 269,555 grains. For the cross- 
sections of three-dimensional microstructures, we took a 
series of cross-sections parallel to the faces of the cubic 
simulation cell. These were spaced roughly five grain di¬ 
ameters apart to reduce correlations between neighbor- 
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ing cross-sections, giving statistics from a total of roughly 
100,000 grains. Though there appears to be no consensus 
in the literature as to the precise fraction of the initial 
grains remaining when the steady-state is reached^, we 
believe our conditions are conservative. 

We report three types of descriptions of steady-state 
microstructures. The first are geometric and topolog¬ 
ical features of individual grains, e.g., distributions of 
grain sizes and number of neighbors. We refer to these as 
point quantities to emphasize that they measure features 
associated with individual grains, as opposed to collec¬ 
tive arrangements of multiple grains. The second type of 
description are correlations of point quantities between 
neighboring grains, where distance between neighboring 
grains is measured between centers of mass using the 
standard Euclidean distance. Third, we investigate cor¬ 
relations of point quantities between neighboring grains, 
where distance between grains is measured using a topo¬ 
logical quantity that we call bond distance; we provide 
evidence that this measurement of proximity is prefer¬ 
able to the more commonly-used notions of nearest and 
next nearest neighbors. Finally, given substantial evi¬ 
dence that the steady-state grain-growth microstructure 
is well-defined, we argue that the steady-state is a more 
reasonable proxy for real polycrystalline materials than 
the frequently-used Poisson Voronoi construction. 

Where possible, we compare analogous data from two- 
dimensional simulations, three-dimensional simulations 
and cross-sections of three-dimensional simulations, de¬ 
noted 2D, 3D, and 3DX, respectively. Error bars in all 
plots indicate one standard error from the mean, though 
any systematic error is neglected. 


A. Point quantities 

The most frequently reported statistics describing 
steady-state grain growth microstructures are of features 
of individual grains. In this section, we consider distri¬ 
butions of point quantities including the number of faces, 
surface area, volume, edge lengths, perimeters, mean 
width, a measure of compactness, and others. A discus¬ 
sion of Lewis’s lav® and Feltham’s lav®, correlations 
between the number of edges and the size of a grain in 
two dimensions, and the generalization of these to the 
three-dimensional structures, is included as well. 


1. Edges and faces 

The von Neumann-Mullins relation indicates that the 
number of edges of a grain in the 2D system completely 
determines the rate of change of a grain’s area, suggest¬ 
ing that the distribution of the number of edges is impor¬ 
tant to the evolution of the micros tructn re as a whole. 
This distribution is reported in Fig. |4(a)| for the 2D and 
3DX systems. Several analogous quantities may be iden¬ 
tified in three dimensions, including the number of edges 



(a) 



(b) 


FIG. 4. (a) Edges per grain in 2D and 3DX, and of faces in 3D. 
Splines provide a guide for the eye; error bars are smaller than 
the markers. Experimental results from a 3DX aluminum 
sampl^^ are provided for comparison, (b) Faces per grain in 
3D; error bars are smaller than the markers. Experimental 
results from a titanium allo}® and a numerical predictior® 
are provided for comparison. 


per f ace or faces per grain, reported in Figs. |4(a)| and 
|4(b) One consequence of Euler’s TheorenP^ is that the 
average number of edges per grain in the 2D and 3DX 
systems must be precisely six. Simulation data shows 
6.000 ±0.001 edges per face in 2D, with a standard devi¬ 
ation 1.273 ± 0.001; analogous data shows 6.000 ± 0.006 
edges per face in 3DX, and a standard deviation of 
1.813±0.006. In 3D, the corresponding constraints on the 
average number of edges per fac e and the average num¬ 
ber of faces per grain is weakei^^I^^. Simulation data 
indicates 5.1283 ± 0.0006 faces per grain in 3D, and a 
standard deviation of 1.1292 ± 0.0006. In 3D, the aver¬ 
age number of faces per grain {F) determines the aver¬ 
age number of edges per face {E) through the relation 
(F) = 12/(6 - (F)). Since (F) = 5.1283 ± 0.0006, we 
have that (F) ^ 13.766 ±0.009. A precise justification of 
these average values in 3D is not known, though several 
studies predict the smaller value of 13.39?^^^^ for the 
average number of faces and one predicts 13.564^. Like- 
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wise, there is no kno wn an alytic function that gives the 
distribut ions in Figs. |4(a)] and (b), despite some recent 
attempt j^^ l ^^ l for the 3D case. Nevertheless, the strik¬ 
ing similarity in the shapes of the distributions for edges 
per grain in 2D and for faces per grain in 3D suggests 
the existence of a shared fundamental cause. Meanwhile, 
the difference in the distributions of edges per grain in 
the 2D and 3DX systems, along with considerable other 
evidence in the literature, makes clear the difficulty of di¬ 
rectly comparing two-dimensional simulations with cross- 
sections of experimental samples. This is the case despite 
suggestive results reported by Anderson et al.^ where 
these distributions nearly overlap. 

Our results for the 2D system compare well with the 
distributio n of edg es per grain as given by other recent 
simnlation d^^ l ^^ l ^ ^ and provide further evidence that the 
distribution peaks at six sides rather than five. This is in 
contradiction to early simulatiorP^ and experimentaP^ 
results. Simulation results for the 3DX case occur less 
frequently in the literature due to the increased compu¬ 
tational requirements, though our distrib ution appears 
to be consistent with more recent studiej^^ ^. Cu ri- 
ously, the corresponding experimental resultP^^^^^^ for 
cross-sections of three-dimensional samples consistently 
give distributions with more five-sided faces than six- 
sided faces. Since the apparent number of edges in cross- 
section depends sensitively on the grain geometry, this 
discrepancy may depend on how a sample was prepared, 
and may not be generically representative of steady-state 
grain growth microstructure. 

For the 3D system, the distribution of edges per face 
in our simulations is consistent with that reported for 
grains in an aluminum allo^l^aswell as with the results 
of several other simulation^^^EHHl^ Finally, the distri¬ 
bution of faces per grain has been a subject of frequent 
study, as it was believed that a grain’s number of faces 
guides its growth in 3D in the way in which a grain’s 
number of edges guides it in 2D. Our findings agree 
with those of simu lations implemented using a variety of 
approache d^^ l ^^ l ^^ l, as well as with the statistics derived 
from a population of thousands of grainP^ obtained from 
a /d-titanium alloy by serial sectioning techniques. There 
is, therefore, ample support for our findings. One of the 
more recent theoretical predictionP^ generally performs 
well, but shows statistically significant deviation from our 
data - it predicts more grains with few faces and fewer 
grains with many faces than we observe. 

Regardless of the similarities between Figs. |4(a)| and 
(b), we note that there is an important difference in what 
they measure. In two dimensions, the combinatorial type 
of a grain is completely described by its number of faces. 
That is, any grain with n edges has exactly the same 
combinatorial structure as any other grain with n edges. 
However, three-dimensional grains are much more com¬ 
plicated, and many distinct polyhedra share the same 
number of face^^. For example, grains with 8 faces 
can have one of 14 distinct types, while those with 9 
faces can have one of 50. As the number of faces in¬ 


creases, the number of distinct combinatorial types in¬ 
creases faster than exponentially. Therefore, the infor¬ 
mation indicated in Fig. |4(b)] is much less descriptive of 
a three-dimensional system than Fig. |4(a)| is of a two- 
dimensional one. 


2. Areas and volumes 


One of the more accessible and physically relevant de¬ 
scriptions of a microstructure is the distribution of grain 
sizes, though that can be measured in different ways. 
The most common measurement of size is the effective 
grain radius, i.e., the square root of the area of a two- 
dimensional grain or the cube root of the volume of a 
three-dimensional grain. Figure [5(a)| gives the distribu¬ 
tions of the effective grain radii for the 2D and 3DX sys¬ 
tems, as well as of the effective face radii for the 3D 
system; the prediction of HillerfP^ for the distribution 
of grain radii in the 2D system is plotted as well for 
comparison. While Hillert’s prediction is still often re¬ 
garded in the literature as a relia ble re ference, it does 
not agree with our results. Figure [5(b)| gives the corre¬ 
sponding distribution for the effective grain radii in the 
3D simulation and from thousands of grains in an experi¬ 
mental /3-titanium allo}!^. The Hillerij^, WeibulP^, and 
RioJ^ distributions are provided for comparison, where 
the Weibull distribution is given by 


f{r) 


kr^ 1 


exp 



( 3 ) 


with fitting parameters k = 2.762 ±0.001 and A = 
1.1570 ± 0.0003 as determined by least squares. The co¬ 
efficients of determination - 0.767, 0.991 and 0.998, 
respectively - indicate that the Rios distribution follows 
our simulations most closely. However, the Weibull dis¬ 
tribution performs only slightly less well and has the ad¬ 
vantage of a simple algebraic expression. _ 

Several features of the distributions in Fig. |5(a)] deserve 
further mention. One is that the distributions for the 
3DX and 3D systems do not pass through the origin, 
meaning that there is a finite probability of finding a 
grain or a face with an arbitrarily small effective radius. 
The other is that the distribution for the effective radii 
of the 2D system shows more than one peak; the most 
pronounced being to either side of an effective radius of 
one. A shoulder occurs just below an effective radius 
of 0.5 as well, though this is more subtle. Given the 
relatively small size of our errors, these features cannot be 
attributed to noise, and another explanation is required. 

We first discuss the observation of finite probability at 
zero size. This discussion is simplified by changing the 
independent variables of the probability density functions 
in Fig.[^ Specifically, a proba bility density function /(r) 
of the radius r = in Fig. 5(a) is converted into the 


corresponding probability density of the area a 

^ 7 - 1/2 

9{a) = 


( 4 ) 
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FIG. 5. (a) Effective grain radii for 2D and 3DX, and ef¬ 
fective fac e rad ii for 3D. (b) Effective gra in radii for 3D. 
PredictionJ^^^ and experimental results^ are provided for 
comparison. Effective radii given in units of average effective 
radius; error bars are smaller than the marker size. 


and is plotted in Fig. 6(a) From this vantage, the pos¬ 


itive probability of finding a grain or a face with an ar¬ 
bitrarily small effective radius in the 3DX or 3D systems 
is simply the result of the probability density function of 
areas in these systems decaying as near the origin. 

The explanation for the peaks in the distribution of the 
effective radii for the 2D system, and the corresponding 
ones in the probability density function of the areas, is 
more involved. For this purpose, we separate the grains of 
the 2D system into classes with a fixed number of edges. 
The probability density functions for the po pulat ions of 
grains in the separate classes appear in Fig. |6(b)[ where 
they may be observed to be sm ooth functions. Adding 
the sep arate curves in Fig. |6(b)| to obtain the distribution 
in Fig. |6(a)| results in peaks precisely where the distribu¬ 
tion for grains with four sides intersects that for grains 
with five, the distribution for grains with five sides in¬ 
tersects that for grains with six, and the distribution for 
grains with six sides intersects that for grains with seven. 
Presumably other peaks occur at the remaining intersec¬ 
tions of these distributions, but the number of grains 


FIG. 6. (a) Grain areas in 2D and 3DX, and face areas in 3D. 
(b) Probability density functions for grain areas in 2D with a 
fixed number of edges. Areas are given in normalized units; 
error bars are smaller than the marker size. 


involved may be small enough that the peak is not easily 
distinguished from the background. This phen omenon 
has been recognized in the literature before^^l^, though 
the apparent importance of the combinatorial structure 
of the grains to the system as a whole does not seem to 
be widely appreciated. This may be due to the relative 
scarcity of sufficiently precise data to adequately resolve 
the peaks, leading to several sug gestio ns that the areas 

3rP 


from Fig. |5(a)| to 
_ ' _ the prob abilit y density function /(r) for the 
3D system in Fig. |5(b)| is converted into a function of 
volume V 


follow an exponential distributio: 

Following the transformation 
Fig. |6(a) 


5(^') = —(5) 

and is plotted in Fig. The most notable feature of the 
resulting probability density function is the resemblance 
to an exponential distri bution. This striking result has 
been reported previousl}^^^^^^^. An exponential distri¬ 
bution of volumes implies a Weibull distribution for ef¬ 
fective grain radii with k = 3 and A = 1 in Equation 
found by inverting Equation]^ The reason why a (near) 
exponential distribution is observed remains unknown. 









































6 



Volume 

FIG. 7. Grain volumes, given in units of average grain volume; 
error bars are smaller than the marker size. 


A second notable feature of the grain volume distribu¬ 
tion is the absence of multiple peaks analogous to those 
for the area distribution of grains in the 2D system. One 
possible explanation is that the separate distributions 
corresponding to those in Fig. 6(b) more nearly overlap!^, 


making the transitions less abrupt. Alternatively, the 
MacPherson-Srolovitz relation indicates that the combi¬ 
natorial type of a grain in a 3D system does not directly 
govern its volume evolution, suggesting that a division 
of the grains into distinct classes based on combinatorial 
type may not be especially meaningful. 

We note that many grain size distributions h ave been 
rep orted (e .g., see recent s imulation data for 2 Dff 31 IIl 2 l^ 
and 3F| 26 | 33 | 67 |^ ^ While these results appear 
consistent between the different studies, the very large 
samples represented in the data presented here bring new 
features to the fore, including the statistically significant 
deviation of the grain volume distribution from a pure 
exponential. 


3. Lewis’s law, Feltham’s law, and generalizations 

While the importance of combinatorially distinct pop¬ 
ulations of grains to the peaks in the grain area distri¬ 
bution function in 2D is not widely recognized, the im¬ 
portance of combinatorially distinct classes has been con¬ 
sidered in other contexts. When investigating a section 
of cucumber epidermis, LewiJ^ noticed that cell areas 
appeared to be proportional to their number of sides n, 
for values of n ranging from about four to seven. This 
relationship is often referred to as Lewis’ law, and plays 
an important part in the application of the maximum 
entropy form alism to the structural properties of cel¬ 
lular network^^SMIlH] widely believed that 

Lewis’ law is a general result, holding over a wide range 
of n (neglecting Lewis’ original qualification), the areas 
of the grains in the 2D system do not appea r to follow 
Lewis’ law over an appreciable range; see Fig. |8(a)] Sim¬ 
ilarly, the same figure indicates that Lewis’ law neither 




FIG. 8. (a) Average area of grains (2D and 3DX) or faces 
(3D) with a given number of edges, (b) Average perimeter of 
grains (2D and 3DX) or faces (3D) with a given number of 
edges. Areas and perimeters given in normalized units; some 
error bars are smaller than the marker size. 


applies widely to grain areas in 3DX nor to face areas 
in the 3D s yst em. This conclusion is supported by other 
literatur^^^^, as well. 


The situation is slig htly b etter for the extension of this 
relation to 3D in Fig. |9(c) , where the average grain vol¬ 
ume appears to increase linearly with the number of faces 
n for n > 20. This is consistent with the literatur^^, 
though Fig. |4(b)| indicates that this extension of Lewis’ 
law would hold only for a relatively small fraction of the 
grains in the system. 


FelthanJ^ provided experimental data on annealed tin 
in support of the claim that the average linear measure of 
a two-dimensional grain increases linearly with its num¬ 
ber of edge s; this is referred to as Feltham’s lawP^^^. 
Figure [8(1^ indicates how perimeters of grains in the 2D 
and 3DX systems, and of faces in the 3D system, increase 
with number of edges. Of course, a single system cannot 
simultaneously follow Lewis’ law and Feltham’s law, and 
the respective accuracy of the two is a question that has 
received some attentiorP^l. With regard to our simula¬ 
tions though, Feltham ’s law does not appear to hold for 
any of the cases in Fig. |8(b)] Nor does Feltham’s law hold 
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(a) (b) (c) 


FIG. 9. (a) Average effective radius, (b) surface area and (c) volume of grains with a given number of faces for 3D. Using least- 
squares fitting, the average volume of grains with n faces for n > 20 follows the linear relation a • n-\-h with a = 0.260 zb 0.001 
and 5 = 2.97 zb 0.03. Effective radius, surface area and volume given in normalized units; some error bars are smaller than the 
marker size. 


for the analogous plot of th e effective radii of grains in 
the 3D system in Fig. |9(a)[ Our statistically significant 
conclusion for an accurate normal grain growth si mula- 
tion is suppo rted by sever al other simulations in 
and 3EI ^^U^ I, but not by alP^^^^^^. 


Finally, Fig. 9(b) [ considers a relationship that lies be¬ 
tween Lewis’ law and Feltham’s law, namely, the depen¬ 
dence of the average surface area of a grain in 3D on its 
number of faces. There is a portion of the plot that could 
be considered linear, though the noticeable deviations for 
low and high values of n discourages any proposal that 
this is a fundamental feature of the system. 


4- Edge lengths, perimeters, and surfaee areas 

Having studied the most common measures of grain 
size in Section III A2[ we now consider th e size of fea¬ 
tures on the boundaries of grains. Figure |lQ(a)| shows 
the distribution of edge lengths for the 2D, 3DX a nd 3D 
systems. As noted previously with regard to Fig. |5(a)[ 
the structure of a genuinely two-dimensional structure 
is noticeably distinct from a cross-section of a three- 
dimensional one. A second striking feature is that there 
is a finite probability of finding an arbitrarily short edge 
in all three systems. However, since the 2D and 3D evolu¬ 
tion equations do not reference the features of individual 
edges, there is no reason to question the existence of ar¬ 
bitrarily short edges in normal grain growth situations. 
Although this distribution appears infrequently in the lit¬ 
erature, ThomaP^ performed a Monte Carlo simulation 
of a 3D system and reported an edge length distribution 
substantially different from ours. We suspect that this is 
attributable to the difficulty of measuring edge lengths 
for voxelized grains, particularly for short edges. 

The absence of the grain perimeter in the evolution 
equations may be used to explain the finite probability 
of an arbitrarily small perimeter for a grain in the 3DX 


system and f or a fa ce in the 3D system as well, as indi¬ 
cated in Fig. |lQ(b)| This does not apply for the pe rime- 
ters of grains in the 2D system though, since Fig. |8(b)| 
indicates that a grain with a small perimeter will gener¬ 
ally have few edges, and therefore becomes smaller (to 
the point o f vanishing). The most conspicu ous fe ature of 
Fig. lQ(b)| though is the resemblance to Fig. |5(a) Indeed, 
they should be related by a rescaling of the independent 
variable, since both give a probability density function of 
a linear measure of the size of a grain. This implies that 
the peaks in the distribution for the 2D system result 
from the separation of the grains into combinatorially 
distinct populations. 

The remaining boundary element is unique to three- 
dimensional grains. Figure shows the distribution 
of surface areas of grains for the 3D system; similar 
distributi ons o f grain surface areas have been reported 
previousl}0^^^. Although the bin for smallest grain sur¬ 
face area shows a substantial probability density, the 
trend in this distribution does not necessarily imply that 
the probability density at zero grain surface area is finite. 


5. Isoperimetrie ratios 

Another way to characterize individual grains is by the 
deviation of their shape from a circle or sphere in two or 
three dimensions, respectively. For a grain in two dimen¬ 
sions or a face in three dimensions, the isoperimetrie ratio 
is defined as the ratio of the area A of a grain to the area 
of a circle with the same perimeter P, or 47rA/P^. Fig- 
ure |12(a)' shows this quantity for grains in 2D and 3DX, 
and for grain faces in 3D. Similarly, for a grain in three 
dimensions the isoperimetrie ratio is defined as the ratio 
of the square of the volume U of a grain to the square 
of the volume of a sphe re havi ng the same surface area 
S, or SGttV^/S^. Figure [r2(b)| shows the distribution of 
this measure for grains from the 3D system. 
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(a) 



(b) 


FIG. 10. (a) Edge lengths in 2D, 3DX and 3D. (b) Grain 
perimeters in 2D and 3DX and face perimeters in 3D. Edge 
lengths and perimeters given in normalized units; error bars 
are smaller than the marker size. 



Grain surface area 

EIG. 11. Grain surface areas, in normalized units; error bars 
are smaller than the marker size. 


The isoperimetric inequality states that the isoperi- 
metric ratio is bounded above by one, with equality 
only for a sphere. Clearly, all of the data peak close to 
unity, indicating that grains in normal grain growth mi¬ 
crostructures are quite close to spheres This is consistent 
with the microstructures seen in Fig. fn. The average 
isoperimetric ratios of grains in the 2D and 3DX sys- 



47rA/P^ 


(a) 



367rVVs^ 

(b) 

FIG. 12. (a) Isoperimetric ratios for grains in 2D and 3DX, 
and for faces in 3D. (b) Isoperimetric ratios for grains in 3D. 
Error bars are smaller than the marker size. 


terns, and of faces in the 3D system, are 0.85628T0.00004, 
0.7817±0.0004 and 0.76803±0.00005, respectively. That 
grains of the 2D system are considerably more circular 
than grains of the 3D and 3DX systems should be ex¬ 
pected from energetic considerations, since the energy 
of the 2D system is proportional to the total boundary 
length. Evolution along the steepest descent in energy 
therefore drives the grains of the 2D system to be more 
circular, subject to the constraint that they remain space 
filling. For comparison, we note that the isoperimetric 
ratio for (space-filling) regular hexagons is ~ 0.907. 

Meanwhile, the average isoperimetric ratio of grains in 
3D is 0.6571 ± 0.0002. This should not be directly com¬ 
pared with the corresponding quantity in two dimensions 
though, since the isoperimetric ratio in two and three di¬ 
mensions is the quotient of quantities with dimensions 
of the second and sixth powers of length, respectively. 
Accounting for this difference by taking the cube root of 
0.6571 gives 0.8694, implying that that grains in three 
dimensions are slightly more spherical than those in two 
dimensions are circular, as may be expected from the in¬ 
crease in geometrical degrees of freedom with dimension. 
For comparison, we note that the isoperimetric ratio for 
(space-filling) truncated octahedra is ~ 0.753. 
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While the energetic argument, above, applies to nor¬ 
mal grain growth structures, the grains of other cellu¬ 
lar microstructures are often substantially less spherical. 
For instance, the cells in Voronoi tessellations of Pois¬ 
son distributed points have average isoperimetric ratios 
of 0.7281 ±0.0001 and 0.53320 ±0.00005 in two and three 
dimensions, respectively. 


6. Mean width and related quantities 


One measure of particular importance in view of the 
MacPherson-Srolovitz relation is the mean width of 
grains in three dimensions. The probability distribution 
for the mean width is plotted in Fig. |13(a) Since the 
importance of the mean width for normal grain growth 
has only recently been appreciateci^, no experimental 
studies and relatively few simulat ions (e .g., sed^ report 
this quantity. Inspection of Fig. |13(a)| reveals tha t this 
distribution nearly coincides with the plot in Fig. |5(b)| 
While the coincidence of these probability density func¬ 
tions does not necessarily imply that the mean width of 
a particular grain scales with the cubed root of its vol¬ 
ume, this observation may help to formulate a closed set 
of equations involving the MacPherson-Srolovitz relation 
{Eq.§, much in the same way that the Hillert theory of 
grain growtlP^ is based upon the von Neumann-Mullins 
relation. 

The sum of the bounding triple edge lengths ap¬ 
pears together with the mean width in the MacPherson- 
Srolovitz relation for grain growth in 3D. To explore 
the interaction of the mean width £ of a grain and the 
sum of the bounding triple edge lengths M. in determin¬ 
ing the growth rate, we consider the averages of these 
quantit ies ove r all grains with a fixed number of faces 
in Fig. |13(b)[ When scaled appropriately, these quan¬ 
tities may be summed to determine the average growth 
rate of grains with a given number of faces (plotted in 
blue). An analysis of regular polyhedrsP^ gives a predic¬ 
tion for the mean width of grains as a function of their 
number of faces, though the difference of the prediction 
from our results reveals that grains with small number 
of faces often differ substantially from regular polyhedra. 
Following the same line of reasoning, a predictiorP^ for 
the rate of change of volume of a grain roughly follows 
the blue curve in Fig. |13(b)| for grains with many faces. 
The difference between this prediction from our results 
for grains with few faces is quite significant though, since 
grains with few faces often evolve the most rapidly. 

Notice that the expected mean width and sum of the 
triple edges of a grain increase with its number of faces 
for grains with many faces. For grains with few faces, 
the mean width dominates the MacPherson-Srolovitz 
equation, and the grain shrinks. For grains with many 
faces, the sum of the triple edges dominates the equa¬ 
tion and the grain grows. The crossover occurs for 
grains with approximately 15 faces. Curiously, there does 
not appear to be agreement between any of the simu- 



(a) 



Faces 

(b) 

FIG. 13. (a) Mean widths of grains, in normalized units; error 
bars are smaller than the marker size, (b) Average mean 
width, average sum of triple edge lengths bounding a grain, 
and average sum of those quantities over all grains with a 
given number of faces. 


lations resultj^* ^^ * ^^ * corresponding to the blue curve in 
Fig. |13(l^ The origin of this inconsistency is not known, 
though the frequent disagreement for grains with small 
numbers of faces suggests that the cause may be an in¬ 
consistent treatment of small grains. 


B. Quantities correlated over metric distance 

Many features of a microstructure depend not only on 
properties of its individual grains, but also on their ar¬ 
rangement relative to one another. A variety of functions 
may be used to characterize these geometric features, in¬ 
cluding t he tw o-point correlation function for two-phase 
material j^^ l ^^ l and the radi al pai r correlation function for 
liquids or amorphous solids^^^^. The quantity measured 
in this section is similar to the usual radial pair correla¬ 
tion function, but is applied to the grains of a microstruc¬ 
ture rather than to atoms of a liquid. 

Specifically, let Xi and yj be values of a property 
of interest (e.g., volume, surface area) for the ith and 
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jf'th grains of the microstructure. The indicator func¬ 
tion Xij{p) is equal to one when the jth grain inter¬ 
sects a sphere of radius p constructed around the center 
of mass of the ith grain, and vanishes otherwise. No¬ 
tice that Xijip) is not necessarily symmetric in i and 
j. Finally, let {x) and (y) indicate the weighted aver¬ 
age values of Xi and yi, calculated using the expressions 

Eij Xij {p)xi/ T,ij Xij {p) and Xij {p)yj / Xij (p), 

respectively. The Pearson correlation coefficient between 
the values of this property for a central grain and the 
grains residing at some distance p is given by 


' xy 


(P) = 


T,Xijip)ixi - {x))iyj - (y)) 
ij 


[T,Xij{p)ixi - (x))^] [T>Xij{p){yj - {y})^] 

ij ij 


( 6 ) 

where the sums include all allowed values of i and j. The 
correlation coefficient Vxy is plotted in Fig. for several 
properties as a function of p. 

The correlation between the numbers of edges of a pair 
of grains for 2D and 3DX, and between the numbers of 
faces of a pair of grains for 3D, is given in Fig. 14(a)| 
We normalize distance by the square root of the average 
grain area or cube root of the average grain volume, as 
appropriate. The most conspicuous feature of this plot is 
the strong peak at small distances, caused by the perfect 
correlation of the number of sides of a central grain with 
itself. As the distance from a central grain increases this 
correlation oscillates, becoming negative at a distance of 
about one average grain diameter and positive at about 
two average grain diameters. This correlation rapidly 
weakens, to the point that there is no noticeable correla¬ 
tion between the number of sides of a central grain with 
grains th ree (av erage) grain diameters or more away. 

Figure |14(b)| shows the correlation between perime¬ 
ters of neighboring grains in 2D and 3DX, and between 
surf ace are as of neighboring grains in 3D. Analogously, 
Fig. |14(c)| gives the correlation between areas of neigh¬ 
boring grains in 2D and 3DX, and between volumes of 
neighboring grains in 3D. Despite the general consistency 
of the behavior described above throughout the plots 


in Fig. 14, a comparison of the correlations in the 2D, 
3DX and 3D systems reveals several significant features. 
For example, the dimensionality of the system appears 
to strongly influence the correlations, to the point that 
the curves for the 2D and 3DX systems nearly overlap. 
This is particularly striking given the differences in the 
statistical features of the 2D and 3DX systems reported 
in Section |mg and the differences in the origin of the 
structures. Meanwhile, the correlations for the 3D sys¬ 
tem are consistently weaker than for the two-dimensional 
systems. This is consistent with the general principle 
that more degrees of freedom are accessible to systems 
in higher dimensions. Finally, the shoulder appearing 
in the plots of the 2D and 3DX systems at a distance of 
about half a grain diameter does not appear to have been 
noticed previously in the literature. 

Some features of these figures may be explained by 


separating the contributions into a part from the central 
grain, a part from th e near est neighbors, and so forth. 
With reference to Fig. |14(a) , at a distance of roughly one 
average grain diameter, the number of sides of a central 
grain is negatively correlated with that of its surrounding 
grains. Assuming that many of the grains intersected by 
a sphere of radius one average grain diameter share a side 
with the central grain, this negative correlation indicates 
a property of the nearest neighbors. Specifica lly, th is 
provides some support for the Aboav-Weaire lavP^^^ or 
the empirical observation that grains with many sides are 
generally surrounded by grains with few sides, and vice 
versa. 

A small positive correlation in the number of sides is 
observed at a distance of two average grain diameters. 
Following the above reasoning, one expects that a sphere 
of this radius should generally intersect the second near¬ 
est neighbors, indicating that the average number of sides 
of the second nearest neighbors is weakly positively cor¬ 
related with the central grain. This may be explained as 
an iterated effect of the Aboav-Weaire law. That is, the 
nearest neighbors of a grain with many sides will gener¬ 
ally have few sides, and the nearest neighbors of any one 
of these nearest neighbors with few sides will generally 
have many sides. Since the set of nearest neighbors of 
the nearest neighbors shares many grains with the set of 
second nearest neighbors of the central grain, the Aboav- 
Weaire law implies that there should be a weak positive 
correlation in the number of sides of a central grain with 
the grains at about two average grain diameters. 

As stated though, the Aboav-Weaire law relates the 
number of sides of grains sharing a particular topological 
relationship, not necessarily a particular spatial separa¬ 
tion. The assumption that the majority of grains at a dis¬ 
tance of one average grain diameter are nearest neighbors 
of the central grain, and that many of the grains at a dis¬ 
tance of two average grain diameters are second nearest 
neighbors of the central grain, appears to be reasonable. 
Nevertheless, the use of a spatial separation introduces 
a complication that could be avoided by considering the 
topological relationship of the grain directly. 


C. Quantities correlated over topological distance 


An interpretation of the correlation functions reported 


in Section |IIIB| relies on the separation of contributions 
from the central grain, from the grains with which it 
shares a side, and from grains further away. This separa¬ 
tion suggests that the features of the structure be studied 
directly for pairs of grains with a given topological rela¬ 
tionship, rather than indirectly for pairs of grains with a 
given spatial separation. Several relevant notions of topo¬ 
logical distance appear in the l iterat ure, including shell 
distanc^^^l®! and bond distanc#^^. The bond distance 
between two grains is the minimum number of edges in a 
path between their boundaries. As a reference. Fig. [^in¬ 
dicates the bond distances of neighboring grains around a 






















13 





FIG. 14. (a) Correlation between the number of edges (2D, 3DX) or faces (3D) of a central grain with those of surrounding 
grains, (b) Correlation between the perimeter (2D, 3DX) or surface area (3D) of a central grain with those of surrounding 
grains, (c) Correlation between the area (2D, 3DX) or volume (3D) of a central grain with those of surrounding grains. 



FIG. 15. Bond distance calculated for neighboring grains in a 
two-dimensional microstructure. The central grain is labeled 
C, and numbers indicate bond distances to other grains. 


central grain in a simulated two-dimensional microstruc¬ 
ture. The shell distance between two grains is the min¬ 
imum number of grain boundaries that must be crossed 
to go from the interior of one to the interior of the other. 

As in Section III B, let Xi and yj be values of a property 
of interest (e.g., the number of sides of a grain) of the ith 
and jth grains. The indicator function Xij{d) is equal 
to one when the jth grain is at bond distance d from 
the ith grain, and vanishes otherwise. The joint proba¬ 
bility f{x^y\d) is calculated from all pairs of Xi and yj 
for which Xij{d) does not vanish. The probability f{x\d) 
is calculated by summing f{x^y\d) over all possible val¬ 
ues of y. The conditional probability for y given x is 
f{y\x^d) = In the following figures, we give the 

mean and standard error of the mean of this conditional 
probability distribution as functions of the value x of the 
central grain and for d = 0, 1 and 2. 

We now consider the number of edges or faces in two 
or three dimensions, respectively. Earlier studies investi¬ 


gated the relationship between the number of sides of a 
central grain to the average number of sides of neighbor¬ 
ing grains (i.e., bond distance 0). This stems from the 
experimental observation by Aboa^^^ on a cross-section 
of polycrystalline MgO that a grain with many sides is 
typically surrounded by grains with few sides, and vice 
versa. Th is equ ation, frequently known as the Aboav- 
Weaire law^^l^, attempts to capture this relationship by 
expressing the average number of sides of a neighbor rrin 
as a function of the number of sides of a central grain 
with n sides: 


= {n) - 


{n)c 


M2 


( 7 ) 


where (n) is the average number of sides per grain in 
the system, /i 2 is the variance of the distribution of the 
number of sides , and a is a fitting parameter depending 
on the struct nr 

The correlation between the number of sides of grains 
and the Aboav-Weaire law appear in Fig. For bond 
distance 0, the fitting parameters a are determined by 
least squares to be a = 1.159 ± 0.004, 0.943 ± 0.007 
and 0.074 ± 0.007 for the 2D, 3DX and 3D systems, re¬ 
spectively. While the resulting predictions agree with 
bond distance 0 averages reasonably well, there is sub¬ 
stantial deviation for central grains with small number 
of sides, especially in 2D and 3DX. Indeed, these devia¬ 
tions should be expected when applying Eq. since the 
usual Aboav-Weaire law neglects several terms that are 
present in a more complete analysiJ^. 

Neighboring grains at bond distance 1 generally dis¬ 
play a behavior that is opposite to that for bond distance 
0, i.e., they tend to have few sides when the central grain 
has few and more sides when the central grain has many. 
By contrast with the grains at two average grain diam¬ 
eters in Section IIIB| though, this may not be explained 
as the iterated effect of the Aboav-Weaire law. Consider 
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Edges of center grain 
(a) 


Edges of center grain 

(b) 


Eaces of center grain 

(c) 


EIG. 16. Average number of edges or faces of neighboring grains as a function of the number of edges or faces of a central grain 
for bond distances 0, 1, and 2; where not indicated, error bars are smaller than the marker size. Predictions of Aboave-Weaire 
for d = 0 are shown as black dashed curves, and are obtained using 112 from Eig. [^and a fitted to the data. 



EIG. 17. All events that change the set of grains at bond 
distance 1 in two dimensions. Numbers indicate the bond 
distance of the grain from a central grain C; +/— signs show 
whether the grain containing the symbol gains or loses an 
edge during the event occurring in the direction of the arrow. 


that both grains at bond distances 1 and 2 are nearest 
neighbors of grains at bond distance 0 in Fig. and 
that the average number of sides of grains at bond dis¬ 
tance 2 is almost independent of the number of sides of 
the central grain in Fig. If the behavior of the grains 
at bond distance 1 could be explained as an iterated ef¬ 
fect of the Aboav-Weaire law, then the same behavior 
should be observed for grains at a bond distance 2. This 
is contrary to the simulation results and some other ex¬ 
planation is required. 

The analysis of the average number of sides of grains 
at bond distance 1 is made more manageable by iden¬ 
tifying all events that can change the set of grains at 
this bond distance. The elementary topol ogical transi- 


tions implemented in our simulations 


give six dis¬ 


tinct events in two dimensions (see Fig. 17) and seven 
distinct events in three dimensions. Although our dis¬ 
cussion focuses on two-dimensional events, we assert that 


the three-dimensional case is analogous. The numbers in 
Fig. [^indicate the current bond distance of a grain from 
a central grain C, while the +/— sign within each grain 
indicates whether the grain has gained or lost an edge as a 
result of the transition. Events (a), (b) and (c) in Fig. 17 
apply an edge-flip to one of the edges in the vicinity of the 
central grain; events (d), (e) and (f) involve the collapse 
of a two or three-sided face. Since the relativ e frac tion of 
such faces in the structure is small (see Fig. |4(a) ), these 
events occur infrequently relative to (a), (b) and (c). 

Statistics shown in Fig. derive from normal grain- 
growth microstructures, in which grains with few sides 
tend to become smaller. In 2D, this is a direct con¬ 
sequence of the von Neumann-Mullins relat ion (E q. [^; 
in 3D this is a result of the trends in Eig. |13(b)| along 
with the MacPherson-Srolovitz relation (Eq. [^. A cen¬ 
tral grain with few sides will, therefore, tend to become 
smaller and shed neighbors that replace the grains at 
bond distance 1, as with event (a) in Eig. 17 Assum¬ 
ing that the neighbors of the central grain most likely to 
be lost are those themselves becoming smaller and losing 
sides, grains at a bond distance 1 around a central grain 
with few sides should have few sides as well. 


A different mechanism is required to explain the obser¬ 
vation that grains at bond distance 1 have slightly more 
than the average number of sides when the central grain 
has many. By reasoning analogous to that for the case of 
a central grain with few sides, a central grain with many 
sides will generally grow and acquire more sides. This 
suggests that event (b) will occur frequently and event 
(a) will occur infrequently. Event (c) should occur at an 
intermediate frequency since this edge-flip does not in¬ 
teract directly with the central grain. Notice that the 
most frequent events (b) and (c) move grains from bond 
distance 2 into the population at bond distance 1. Erom 
Eig. grains at bond distance 2 should have the same 
expected number of sides as the system average. While 
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(a) 



(b) 



Volume of center grain 

(c) 


FIG. 18. Average size of neighboring grains as a function of size of a central grain (in normalized units) and bond distance. 
Spline curves have been added to provide a guide for the eye. 


event (b) moves these grains to bond distance 1 with¬ 
out changing the number of sides, event (c) increases the 
number of sides of the relevant grain by one. The average 
number of sides of a grain at bond distance 1 from a cen¬ 
tral grain with many sides should therefore be expected 
to be slightly higher than the system average. 

We now examine the sizes of grains separated by a 
particular bond distance, as indicated in Fig.[^ A con¬ 
sistent feature of these results is that the average area 
or volume of neighboring grains at bond distance 0 de¬ 
creases with increasing area or volume of a central grain; 
this is similar to the behavior of the average number of 
sides of grains at bond distance 0 as a function of the 
number of sides of a central grain in Fig. This sim¬ 
ilarity should be expected given the strong correlation 
between the number of sides and the size of a grain as 
seen in Figs. and More striking, however, is the ob¬ 
servation that for almost any bond distance and any size 
of central grain, the average size of surrounding grains 
is higher than the corresponding system average (i.e., 
unity). This phenomenon is a consequence of the fact 
that grains with large areas and many sides are neigh¬ 
bors to many grains and are included in the averages of 
Fig. more frequently than grains with small areas and 
few sides. This shifts all of the curves upward; this effect 
increases with increasing connectivity of the network as 
occurs in going from two to three dimensions. 

However, this does not explain why the average size 
of neighboring grains at bond distance 1 is consistently 
higher than the average size of grains at bond distance 
2. For this purpose, consider the grains in the vicinity of 
the central grain C in Fig. Suppose that the number 
of sides of the central grain is n. When the central grain 
is surrounded by grains all with six sides (e.g., see (b)), 
the number of grains at bond distance 1 is n and the 
number of grains at bond distance 2 is n. A five-sided 
grain at bond distance 0, as in (a), reduces the number 
of grains at bond distance 2 by one. Conversely, a grain 




FIG. 19. A central grain C with neighboring grains at various 
bond distances. Illustrations highlight the effect of the num¬ 
ber of sides of grains at bond distance 0 on the total number 
of grains at bond distance 2, as described in the text. 


with seven or more sides at bond distance 0, as in (c), 
increases the number of grains at bond distance 2 by one. 
Now, grains of small area generally have few sides (see 
Fig.[^ and grains with few sides generally have neighbors 
with many sides (shown in Fig. 16). A small central grain 
will therefore have more neighbors at bond distance 2 
than at bond distance 1. By analogous reasoning, a large 
central grain will have fewer neighbors at bond distance 
2 than at bond distance 1. Finally, notice that bond 
distance is a symmetric function. Therefore, a large grain 
will be observed more often at bond distance 1 from a 
central grain than at bond distance 2, while a small grain 
will be observed more often at bond distance 2 from a 
central grain than at bond distance 1. This difference in 
sampling will drive the curves in Fig.[^for bond distance 
1 higher than the curves for bond distance 2. The same 
reasoning pertains to the three-dimensional case. 

The final property of interest considered is the perime¬ 
ter or surface area of a grain; results are presented in 
Fig.[^ Overall, the behavior is similar to that in Fig. [18] 
and analogous explanations for this behavior apply. One 
significant difference, however, is the appearance of dis¬ 
tinct peaks in the average perimeter of grains at bond 
distance 0 from a central grain in the 2D syst em. The 
resemblance of these peaks to the ones in Figs. |6(a)| and 
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(a) (b) (c) 


FIG. 20. Average boundary size of neighboring grains as a function of boundary size of a central grain (in normalized units). 
Some error bars are smaller than the marker size. 


|lQ(b)| suggests that they arise from the same source, i.e., 
from the separate contributions of grains with different 
number of sides. 

The results in this section suggest that analysis of a 
microstructure by shell distance is generally not suffi¬ 
cient to distinguish the arrangement or the characteris¬ 
tics of grains around a central grain, and that a more 
complete description of topological relative position is 
required. While bond distance provides one such alter¬ 
native, it is clear that it remains an incomplete measure 
of topological relative position. 


D. Poisson—Voronoi microstructure 

While we consider the steady-state, isotropic grain- 
growth (curvature flow) microstructure to be the pro¬ 
totype or canonical form of a polycrystalline microstruc¬ 
ture, other relatively simples types of microstructures are 
also discussed in the literature. In general, different evo¬ 
lution or microstructure formation laws yield microstruc¬ 
tures with very different statistical (geometric and topo¬ 
logical) properties. Probably the most widely discussed 
“other” microstructure is the Poisson-Voronoi (PV) con¬ 
struction. In a materials science context, the PV mi¬ 
crostructure is often described as arising from the simul¬ 
taneous nucleation of grains at random positions in space, 
and the isotropic growth of these grains such that the ve¬ 
locity of the growth front is constant in time. When two 
growth fronts impinge, they are replaced by stationary 
grain boundaries, resulting in a static PV microstruc¬ 
ture. This microstructure may apply in the ex treme 
limits of phase transformations or recrystallizatiorP^^^. 
Other si milar nucleation and growth laws have also been 
propose d^^P^ l 

To illustrate several simple examples of the differences 
between the PV and the steady-state, isotropic grain- 
growth microstructures, we compare the distribution of 


faces per gra in and grain volume in the two microstruc¬ 
tures in Fig. |21(a)| Not only is the average number of 
faces per grain noticeably smaller in the grain-growth mi¬ 
crostructure, the distribution is also considerably wider. 
The distribution of grain volumes in the two systems is 
also significantly different. Whereas cell volumes follow 
a near-exponential distribution in the grain growth mi¬ 
crostructure, small grains in PV microstructures are van¬ 
ishingly rare and the grain volume distribution is peaked 
at finite grain size. Much more extensive sets of PV 
statistics to compare with the other measurements in this 
paper can be found elsewher^^. 


IV. DISCUSSION AND CONCLUSION 

The existence of a statistically steady-state grain- 
growth microstructure has been commonly conjectured 
in studies of grain growth. To achieve such a statisti¬ 
cal steady-state, the microstructure must undergo grain 
growth for sufficient time for the system to “forget” its 
initial state. In simulations, the initial state is typically 
a Voronoi tessellation of a random distribution of points, 
while in experiment it is inherited from prior process¬ 
ing. Substantial evolution requires growth growth sim¬ 
ulations to begin from a large number of initial grains 
such that when statistics are measured, most of the ini¬ 
tial grains have been consumed. In the present study, we 
report statistical analysis of some of the largest two- and 
three-dimensional normal grain growth simulation data 
sets produced to date. These simulations are accurate in 
that they are consistent wit h the von Neumann-Mullins 
relation in two dimensionJ^^EHl ^^d the MacPherson- 
Srolovitz relation in three dimensionJ^. A detailed de¬ 
scription of the simulation m ethods and their numerical 
errors can be found elsewher d^^ l ^^ l ^^ i 

The simulation results presented here vividly il¬ 
lustrate a number of qualitative differences between 
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(a) 



Volume 


(b) 

FIG. 21. Distribution of (a) faces per grain and (b) normal¬ 
ized grain volumes from grain-growth and Poisson-Voronoi 
microstructures. 


two-dimensional systems and cross-sections of three- 
dimensional ones. The distribution of edges in the cross- 
sections is noticeably wider than in the genuine two- 
dimensional system, as is the distribution of edge lengths. 
Cross-sections thus contain a larger fraction of grains 
with many sides, as well as a higher fraction of very long 
edges. Moreover, grains in the genuine two-dimensional 
system are significantly more circular than grains in 
cross-sections. These three differences can be accounted 
for by considering that two-dimensional systems evolve 
through a steepest descent path in energy, which in two- 
dimensional systems is proportional to the total length of 
the grain boundaries. This inhibits the development of 
very long edges, grains with many edges, and grains with 
low isoperimetric ratios. Because the evolution of cross- 
sections is not governed by the same energy-minimization 
considerations, grains tend to be less circular, and long 
edges and grains with many sides arise more frequently. 

The more prominent feature distinguishing genuine 
two-dimensional systems from cross-sections of three- 
dimensional ones can be observed in the distributions 
of grain areas and perimeters. Two-dimensional sys¬ 
tems exhibit bimodal distributions of these quantities. 


whereas cross-sections exhibit unimodal ones. Although 
these features have been observed previously, insufficient 
data have been available to draw definitive conclusions. 
Given the accuracy in the data reported here, these two 
trends cannot be attributed to noise. They appear to re¬ 
sult from the role that topology plays in the evolution of 
genuine two-dimensional systems, which it does not play 
in the evolution of cross-sections. 

We also observe important qualitative differences be¬ 
tween the 2D and 3D systems. In particular, we note 
that unlike the bimodal distribution of grain areas and 
perimeters in two-dimensional systems, the distribution 
of grain volumes and surfaces areas in three dimensions 
is unimodal. This may similarly be attributed to the 
reduced role of topology in three dimensions: whereas 
the von Neumann-Mullins relationship relates the growth 
rate of a grain to its topology, the analogous MacPherson- 
Srolovitz relation does not. 

Another significant finding of our simulations is the 
near-exponential distribution of grain volumes in three- 
dimensional systems. The probability density function 
of finding a grain with a given normalized volume x is 
accurately described by the simple exponential e~^. The 
reason for this near-exponential distribution of the grain 
volumes is not known and requires further inquiry. 

We also find that while neither Lewis’ law or Feltham’s 
law accurately describe the relationship between the 
number of edges of a grain in two dimensions and its 
size, an similar rule seems to hold in three dimensions. 
That is, the volume of a grain appears to grow linearly 
with its number of faces n for n > 20. It should how¬ 
ever be noted that the fraction of grains with 20 or more 
faces is rather small. We do not currently understand 
the source of this relationship. 

The last important set of results concerns pairwise cor¬ 
relations between grains separated by some metric or 
topological distance. Our results illustrate a qualitative 
similarity between two-dimensional systems and cross- 
sections of three-dimensional ones. The pairwise correla¬ 
tions between grains separated by some metric distance 
appear very similar in both systems despite the substan¬ 
tial difference in their origins. Both of these systems 
exhibit correlations noticeably stronger than those that 
appear in the three-dimensional system. It appears that 
the dimension of the systems plays a larger role in deter¬ 
mining these correlations than the particular method of 
generating them. In all three systems, we find that pair¬ 
wise correlations of sizes and shapes are roughly zero at 
distances greater than three grain diameters. This indi¬ 
cates that although the rules of evolution impose some lo¬ 
cal ordering on a structure, this ordering vanishes rather 
quickly. 

Pairwise correlations between grains at different topo¬ 
logical separations provide additional insight into grain 
growth structures and methods of measuring their prop¬ 
erties. Neighboring grains at bond distance 0 tend to 
have fewer sides when a central grain has more sides and 
vice versa; this is consistent with the classical Aboav- 
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Weaire law. When considering generalizations of the law 
to grains that are not immediate neighbors, we observe 
that data for bond distance 1 is almost always distin¬ 
guished from data for bond distance 2. This suggests 
that the traditional notion of next-nearest neighbors - 
which fails to distinguish between these two distances - 
is inadequate in analyzing grain-growth microstructures. 
In particular, we find that grains at bond distance 1 tend 
to have few sides when the central grain has few sides 
and vice versa] this correlation does not hold for bond 
distance 2. All correlations between a central grain and 
its neighbors vanish for grains at bond distances 3 or 4, 
which illustrates that the ordering generated by grain- 
growth evolution is fairly local. 

We also consider generalizations of the Aboav-Weaire 
law to correlations between metric properties of neighbor¬ 
ing grains. We note that for bond distance 0, the aver¬ 
age area, perimeter, volume, and surface area of a neigh¬ 
boring grain decreases with increasing area, perimeter, 
volume, and surface area of a central grain. This can be 
explained by the strong correlation between the topology 
of a grain and its metric size as investigated earlier in the 
paper. This data also illustrates the role which topology 
plays in the genuine two-dimensional structure, which ex¬ 
hibits multiple peaks in the data for correlations of areas 


and perimeters between neighboring grains; these peaks 
do not appear in the data from the cross-sections. 

Finally, since the steady-state grain-growth mi¬ 
crostructure is reasonably well-defined, a comparison 
with the Poisson-Voronoi microstructure is meaningful. 
The Poisson-Voronoi microstructure is found to be dras¬ 
tically different, suggesting that it is a poor representa¬ 
tion of polycrystalline microstructures that have expe¬ 
rienced extensive grain growth. Nonetheless, the ease 
of constructing PV microstructures has led many re¬ 
searchers to use them to represent polycrystalline materi¬ 
als independent of the physical contexfP^HSSl^ respect¬ 
fully argue that this practice should be discontinued. 
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